Kinetics of a superconductor excited with a femtosecond optical pulse 
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Superconducting state dynamics following excitation of a superconductor with a femtosecond 
optical pulse is studied in terms of a phenomenological Rothwarf and Taylor model. Analytical 
solutions for various limiting cases are obtained. The model is found to account for the intensity 
and temperature dependence of both photoinduced quasiparticle density, as well as pair-breaking 
and superconducting state recovery dynamics in conventional as well as cuprate superconductors. 



In recent years numerous studies of non-equilibrium 
carrier dynamics in superconductors (SC) have been per- 
formed utilizing femtosecond real-time techniques [1-18]. 
Research focused on the identification of relaxation pro- 
cesses and direct measurements of the relaxation times. 
One of the open issues at the moment is vifhether cuprates 
are in the so called phonon bottleneck regime as conven- 
tional SCsJ^j, or in the weak bottleneck regime, where 
relaxation is governed by the biparticle recombination 
kinetics 0, The theoretical model that has 

been most commonly used to interpret the dynamics is 
a phenomenological Rothwarf- Taylor (RT) model which 
describes the evolution of quasiparticle (QP) and high 
frequency phonon (HFP) populations via a set of two 
non-linear differential equations 20], which were shown 
recentlyjl^ to follow from the general set of kinetic equa- 
tions for a SC[23]. While the RT model has been known 
for almost 40 years, no rigorous attempt to solve it has 
been made thus far, and neither has a comparison to the 
experimental data been made. 

In this Letter we present a detailed study of the evo- 
lution of the SC state following excitation by ultrashort 
laser pulse using the RT model. We have obtained an- 
alytical solutions of the model in the limit of a strong 
and a weak bottleneck, which are in excellent agreement 
with numerical simulations. The solutions enable com- 
parison of the model to the experimental results. We 
show that RT model can account for most of the exper- 
imental observations in conventional as well as cuprate 
SC, and that both conventional and cuprates SCs are in 
the strong bottleneck regime, where SC state recovery is 
governed by the HFP decay dynamics. 

Rothwarf and Taylor have pointed out that the phonon 
channel should be considered when studying the SC 
relaxation|20]. When two QPs with energies > A, where 
A is the SC gap, rccombine a HFP {uj > 2A) is created. 
Since HFP can subsequently break a Cooper pair creat- 
ing two QPs the SC recovery is governed by the decay 
of the HFP population. The dynamics of QP and HFP 
populations is determined by|2flj: 

dn/dt = /o + r]N - Rn"^ 

dN/dt^ Jo-T]N/2 + Rn^/2--f{N ~ Nt) (1) 
Here n and N are concentrations of QPs and HFPs, re- 



spectively, J] is the probability for pair-breaking by HFP 
absorption, and R the bare QP recombination rate with 
the creation of a HFP. Nt is the concentration of HFP in 
thermal equilibrium at temperature T, and 7 their decay 
rate. Iq and Jo represent the external sources of QPs 
and HFPs, respectively |0]. Physically 7 is governed 
by the fastest of the two processes: anharmonic decay 
of HFP {uj < 2A phonons do not have sufficient energy 
to break Cooper-pairs) Q and the diffusion of HFP into 
the substrate uSL, While the main T-dependence in 7 
appears near Tc;4|, some T-dependence is expected at 
low-T as well. However, as shown in Refs.|22l |2J] this 
T-dependence is very weak in the temperature range of 
interest and we can consider 7 to be T- independent. 

Since rj/R has the dimensionality of concentration we 
introduce dimensionlcss QP and HFP concentrations, 
q = Rn/r] and p = RN/rj, while 9 = rjt and 7 = 7/77 are 
the dimensionlcss time and HFP decay rate, respectively. 
Since a femtosecond optical pulse is usually shorter than 
characteristic timescales of SC dynamics, Iq and Jq could 
be approximated by (^-functions leading to the initial con- 
centrations of QPs qo — Rno/rj and HFP po = RNq/j] 
p^ . As a result Eqs. |^ reduce to 

dq/de =p-q^ (2) 
dp/de = -p/2 + q^/2 - 7(p - pt). (3) 

with the initial conditions p(0) — pq, q{0) = go- Here 
Pt — RNt /rj is the dimensionlcss concentration of HFPs 
in thermal equilibrium. Thermal equilibrium concentra- 
tions of HFPs and QPs (qt) satisfy the detailed balance 
equation pT = q^- 

We investigated various limiting cases. We refer to 
the situation when the photoinduced density is small 
(po — Pt), (90 — Qt) ^ 1 as a weak perturbation regime, 
while the opposite case is referred to as a strong pertur- 
bation. In case 7 <C 1 we have a strong bottleneck (quasi- 
equilibrium between QPs and HFPs is established, and 
SC state recovery is governed entirely by 7) while in the 
opposite case (7 ^ 1) SC is in a weak bottleneck regime. 

It is useful to estimate the values of thermal QP and 
HFP concentrations and compare them with rj/R. If 
D{ll!) = 9h'LL!'^ / LOjj is the phonon density of states with 
h' the number of atoms per unit cell and lu]j the Debye 
energy, it follows that Nt = ^e^Afr exp(-2A/fcBr). On 
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the other hand the QP density per unit cell at tempera- 
ture T is given by "1: 



ut N{0)y/2TTAkBTexp{-A/T). 



(4) 



Here A^(0) is the electronic density of states per unit cell. 
From rj/R = u'^/Nt follows that r]/R 



JV(0)^ 



Re- 



gardless of the value oi -q/R, it follows that the high 
perturbation limit is reached when the photoexcitation 
density is close to the density required for complete de- 
pletion of the SC state. To show this, we estimate the 
temperature T where qf ^ 1 (uf sa "u/R )■ It follows 
that f w A/ln( ^^j/^ ); thus nr is comparable to rj/R 
only in the close vicinity of T^- 

Strong bottleneck regime. Let us consider the so- 
lution of Eas. (|2l3|l for the case when 7^1. In this 
case two distinct regimes can be defined: (I) prebot- 
tleneck dynamics, describing the processes on timescale 
much shorter than 7"^, and (II) the SC state recovery at 
9 > 1/7. 

(I) Prebottleneck dynamics describes the short 
timescale evolution of concentrations p{9) and q{9) pre- 
ceding the relaxation. For such short times we can ne- 
glect the last term in Eq. Q . It leads to the conservation 
law q{9) + 2p{9) ~ qo + 2po and exact solutions for QP 
and HFP populations, respectively, in the pre-bottleneck 
regime are obtained|l2|: 
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1 - Keiip{-9/C) 

^-2 ^-1 



l-Kexp {-9/S) 



(5) 
(6) 



Here = ^l/4 + 4po + 2go and K = ^t^tlKW-^ ' 
times 9 > 1 concentrations q{9) and p{9) reach the quasi 
stationary solution 
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qs = -(v/l + 16po + l 
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(7) 



The prebottleneck dynamics depends on the initial con- 
ditions, and has two distinct regimes characterized by the 
parameter K. The regime < K < 1 corresponds to the 
situation when go > 9s, and q{9) rapidly decreases during 
the formation of the bottleneck - Fig. 1(a). On the other 
hand ~l < K < represents the situation when po > Ps, 
and the QP density increases up to times 6 - Fig. 

1(b). This situation is realized in MgB2 0|. 

(II) Superconducting state recovery dynamics. At large 
times {6 > I/7) right hand side in Eas. H2l3|) is determined 
by the HFP decay term. To describe the recovery we note 
that at large times the difference q^ — p « l{p — Pt) 
while p(9) and q{9) are slowly decreasing. We intro- 
duce a function s{9) = q{9Y — p[9). Cleary s{9) <^p,q 
and is slowly decaying {ds/d9 <C s). By neglecting the 
derivative ds/d9 it follows from the definition of s{9) that 
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^, which simplifies Eas. H2l3() to 



(l + 27 + 4g(0)) 



dm 

d9 



-2^{q{9 f - 4) 



(8) 



After integration, we obtain an analytical solution: 

1 + 27., ,q~qt. 1 -f 27 ,q + qT 

1 m 



^2761 = (2- 



2qj 



<ls - qr 



2qT 



Qs + qr 
(9) 

Combining solutions for prebottleneck and recovery dy- 
namics we obtain an approximate solution for q{9) valid 
over the entire time range 



qi9) ^ qii9) + q2i9) 



(10) 



Here qi{9) and q2{9) are determined by Eqs.© and © 
respectively. To illustrate that Ea. (|10|l describes the so- 
lutions of Eas. (|2l3|l we compare the numerical solution 
with approximation given by Ea. Hl()|l . Fig.l shows excel- 
lent agreement for all the limiting cases. 




FIG. 1: Analytical solutions of RT equations (open symbols) 
compared to numerical solutions (solid lines) in several limit- 
ing cases at T = K. Panels (a) and (b) represent the strong 
bottleneck limit (7 = 0.1) for different initial conditions (see 
text). Analytical solution for ^(6^) is given by Eq.lCni. Panel 
(c) presents the weak bottleneck case (7 = 10). The inset to 
panel (a) presents the inverse of the photoinduced QP den- 
sity as a function of 6 for various temperatures (qr)- At low- 
T {qr — > 0) the time dependence is linear, mimicking the 
bimolecular kinetics, while at higher T it becomes nearly ex- 
ponential (dashed curve presents an exponential decay fit). 
Indeed, similar behavior was recently reported fr§| . 



Eq.0 can be simplified in the low and high-T limits: 

i) Low temperature limit (qx ^1). 

In this case 2 ^ and Eq.lO reduces to: 



q{9) = qr 



^qrjqs " qr) exp {-9/tt) 
{qs + qr) - {qs - qr) exp {-9/tt) 



(11; 
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The relaxation rate = 475^/(1 + 27) is strongly re- 
duced by the additional T-dependent factor qx cc ut- 
The behavior is consistent with the T-dependence of the 
relaxation time r in most of the SCs: cupratej3, 0, 
and conventional 0, The low-T divergence of r 

was origiiially attributed to the bi-particle prebottleneck 
dynamics 9] where r is divergent due to strong reduction 
of the thermal QP density. Here we show that at low-T r 
diverges also in the case of a well established bottleneck. 

In the limit of 0/tt <C 1 Ea. Hll(l further reduces to 
q{9) — gs/(l + showing identical dynamics to 

that of the biparticle recombination - see inset to Fig. 
1(a). If one defines the relaxation rate as the slope 
= \dq{9)/d9\/{qs — qr) at ^ 0, as proposed in 
Ref.0|, the SC relaxation rate is intensity dependent 



2(gs + qx)! 

(1 + 27) ■ 



(12) 



Indeed, the intensity dependent recovery dynamics has 
been recently observed at low-T in BiSC0[l3|, in agree- 



The intensity dependence is observed only if qs^ qr] in 
this case is proportional to the excitation intensity. It 
follows from Ea. H12|l . that the observation of pronounced 
intensity dependence of r is constrained only to the very 
low-T (determined by the A/kTc ratio) and is more likely 
to be observed in cuprates than in conventional SCs. 

overy dyna 
JiSCoQ, 

ment with the solution of RT equations in the strong 
bottleneck regime - inset to Fig. 1(a). 

ii) High temperature limit {qx Jc^ 1). 

In this case 2 ^ ^^'^ Eq.® reduces to: 

q{e) = y^<z| + (<?2-q2)exp i-e/TT) (i3) 

Here the relaxation rate t^^ = 7 is weakly T-dependent 
(due to intrinsic T-dependence of 7), and does not de- 
pend on qs or ps, implying that the relaxation rate is 
intensity independent. The same is true also for t~^. 

Weak bottleneck (7^1). In this case at short 
times 9 « 7"^ ^ 1 the last term in Eq.© is dominant, 
and the solution has the form 



p{9) ^pt + (pq - pt) exp (-76I) 



(14) 



The HEP density reaches its thermodynamic value on 
a time scale determined by 7~^. On such a short time 
scale QP density has not been changed, therefore we can 
substitute p with pT in Eq.ll^J. Taking into account that 
9t ~ Pt obtain the solution: 



gT[go + gT + (go - qr) exp {-~2qT9)] 
qa+qr - (go - qr) exp {-2qT9) 



(15) 



Eas. H14ll5|l present the solution of Eas. (|2l3(l for aU 

timescales. As seen in Fig. 1(c) the difference between 
the analytical and numerical solutions is negligible. Sim- 
ilar to the case of a strong bottleneck the relaxation rate 



determined as the slope of the signal at = is found to 
be intensity dependent if go S> gT with = {qg + qx). 
The main difference between the strong and weak bot- 
tleneck cases is that the absolute value of is in the 
strong bottleneck case reduced by the HEP decay time. 

Photoinduced QP density. Earlier 0] we have dis- 
cussed the T and excitation intensity dependence of the 
photoinduced signal amplitude Q (proportional to the 
photoinduced QP density) assuming the bottleneck con- 
dition and the energy conservation law. These results 
were recently confirmed using a realistic phonon den- 
sity of states[2^. However, solution of Eqs.|2El allows 
unambiguous comparison of Q{T) with the expected T- 
dependence in a SC - without an additional assumption 
of the energy conservation as in Refs. 0, . 

In the strong bottleneck regime initial dynamics leads 
to the stationary densities of QPs and HFPs determined 
by Eqs.Q. The initial concentrations can be written as 
go = gT + and po = pr + Ap, where Ag, Ap are 
photoinduced initial concentrations of QPs and HFPs, 
respectively. Assuming weak perturbation limit Ap, 
Aq <^ 1 the amplitude of the signal is given by: 



Q{T) acqs - qr 



2Ap + Aq 



Ay/1 + 16pT 



(16) 



Taking into account that px — g^ and gT = RriT/r], 
and normalizing Q to its low-T value, A — Ql Qt^qki 
we obtain ut oc A^^ — 1. This relation is expected to be 
general and valid for any SC irrespective on the gap sym- 
metry allowing direct estimation of the T-dependence of 
QP density in thermal equilibrium. 

Importantly, in the weak bottleneck case Q should be 
proportional to the pump intensity, and T-independent 
at low-T since q{9 = 0) = go, while increasing when 
T ^ Tc due to closing of the SC gap. This T-dependence 
has not been observed thus far. 

In Fig. 2 we plot the T-dependences of ht (extracted 
from the raw data via nr oc A^^ — 1) and the SC 
recovery time t for MgB2 and the cuprate supercon- 
ductor Tl2Ba2Cu06+5. The two datasets capture the 
main features observed on other conventional 0, Q| 
or cuprate 0, IE S S superconductors. In both 
gradually decreases upon increasing T towards Tc 
, 0, II H H Elli , while the low-T relaxation time 
increases when T — > 0. The only exception thus far 
was the near optimally doped YBCO, where r is roughly 
constant at low-T 0, 0. The main difference between 
cuprates and conventional SC is a 1-2 orders of magni- 
tude difference in t, which can be understood in terms of 
an order of magnitude difference in the value of the SC 
energy gapfl^. 

While at low-T and weak perturbations both strong 
and weak bottleneck scenarios suggest slowing down of 
the relaxation upon cooling as well as intensity depen- 
dent relaxation, the fact that the amplitude gradually 
decreases upon T-increase suggests that all SC studied 
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thus far are in the strong bottleneck regime (in the weak 
bottleneck scenario Q should increase upon increasing 
T). In order to test the agreement of the model with the 
data, we fit tit with Eq.g)), where BCS form for A(T) 
was used. We find good agreement between the extracted 
values for 2A/Tc from the fit with the values from other 
experimental techniques 0, |2^. 
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using high repetition rate laser systems like in Ref.|5|. 
Finally, we should note that the value of A extracted is 
also somewhat dependent on the T-dependence of A, i.e. 
when fitting we used a BCS form. However, the accu- 
racy of the available experimental data is not sufficient 
enough to speculate about the deviation of the tempera- 
ture dependence from the BCS functional form. 

In conclusion, we studied the evolution of a SC state 
following perturbation with ultrashort optical pulses. We 
derived analytical solutions of RT model for the strong 
and weak bottleneck regimes. Analytical solutions ac- 
count for both pair-breaking dynamics 12] as well as SC 
state recovery, and are in excellent agreement with nu- 
merical results. They enable comparison to experimen- 
tally measured quantities like the SC state recovery time 
and the photoinduced QP density. Comparison with pub- 
lished data suggest that both cuprate and conventional 
SCs follow the RT scenario, both being in the strong bot- 
tleneck regime. 

We would like to acknowledge A.S. Alexandrov and T. 
Mertelj for critically reading the manuscript. 



FIG. 2: Comparison of the experimental data on T- 
dependence of tit and r on a) MgBa [H m and b) 
Tl2Ba2Cu06+5 0] with the RT solution, tit is extracted from 
the raw data (shown in insets) using ht oc A~^ — 1 and fit with 
Eq.||lJ, while r was fit with Ea. H12^ . The extracted values of 
2A/Tc ratio are also shown. 

We fit the T-dependence of the relaxation time r with 
Ea. (|12|l . Assuming 7 to be T- independent and express- 
ing qr with Eq.(|3I), Eq.lO is rewritten to t(T) = !/($ + 
A(AT)i/2exp(-A/T)). Here $ = B{q, - qr), while A 
and B are fitting parameters. Since the T-dependencc of 
qs — qr is measured directly - see Eq. , we are left with 
3 fitting parameters: A, B and A. The T-dependence of 
r is governed mainly by exp(— A/T) term, while A and 
B account for the magnitude of qT and {qs — qr) respec- 
tively. At the intermediate temperatures qs — qr Or, 
and r is governed by the T-dependence of gr showing 
exp(A/T) behavior. At low enough T, however, qs » qr 
and the relaxation time saturates. The values of A ex- 
tracted from relaxation time data are somewhat lower 
than the values extracted from the fit to ut- This can be 
partially attributed to the fact that the T-dependence of 
7 was neglected. More importantly, improper (or miss- 
ing) treatment of the laser heating can also cause an un- 
derestimate of A obtained from the fit. Namely, due to 
laser excitation the equilibrium T of the probed spot can 
be substantially higher than that of the cold finger |2^. 
especially at low-T. Therefore, the T-dependence of r is 
not as steep as if the temperature scale was corrected 
for heating, giving rise to an underestimate of A. Note 
that heating is particularly pronounced in experiments 
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